"Clusterization" and intermittency of temperature fluctuations in turbulent 
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Temperature time traces are obtained in turbulent thermal convection at high Rayleigh numbers. 
Measurements are made in the midplane of the apparatus, near the sidewall but outside the boundary 
layer. A telegraph approximation for temperature traces is generated by setting the fluctuation 
amplitude to 1 or depending on whether or not it exceeds the mean value. Unlike the standard 
diagnostics of intermittency, the telegraph approximation allows one to distinguish the tendency 
of events to cluster ("clusterization") from their large-scale variability in amplitude. A qualitative 
conclusion is that amplitude intermittency might mitigate clusterization effects. 
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I. INTRODUCTION 

We consider turbulent convection in a confined con- 
tainer of circular crosssection and 50 cm diameter. The 
aspect ratio (diameter/height) is unity. The sidewalls are 
insulated and the bottom wall is maintained at a constant 
temperature, which is higher by a small amount AT than 
that of the top wall. The working fluid is cryogenic he- 
lium gas. By controlling the temperature difference be- 
tween the bottom and top walls, as well as the thermo- 
dynamic operating point on the phase plane of the gas, 
the Rayleigh number (Ra) of the flow could be varied 
between 10^ and 10^^. We measure temperature fluctua- 
tions at various Rayleigh numbers towards the upper end 
of this range, in which the convective motion is turbulent. 
Time traces of fluctuations are obtained at a distance of 
4.4 cm from the sidewall on the center plane of the ap- 
paratus. This position is outside of the boundary layer 
region for the Rayleigh numbers considered here. More 
details of the experimental conditions and measurement 
procedure can be found in Ref. 0. 

A significant part of convection, even at the high 
Rayleigh numbers that concern us here, is due to plumes 
01 • We use the term here merely to denote an orga- 
nized activity of convection without implying much about 
their three-dimensional shapes and sizes, or the param- 
eters on which they scale, though a few comments will 
be made momentarily. The primary goal of the paper 
is to learn about the tendency of the plumes to cluster 
together ( "clusterization" ) . 

The upper part of Fig. 1 shows a short segment of tem- 
perature fluctuations at Ra = 1.5 x 10^^. There are four 
large-scale events within this segment (marked by the 
letters A-D), and we imagine them to be the manifesta- 
tion of large-scale plumes. Each event consists of several 
subevents, and there also exist a number of small events 
marked a-i. While it may well be that the subevents 
deserve to be considered separately, we regard them col- 
lectively here. Under these circumstances, it is clear that 
a typical life-time of the large events is of the order of 8 



seconds. Noting from Ref. Q that the mean speed of the 
large-scale circulation ( "mean wind" ) for these conditions 
is about 6 cm/s, a typical length scale of these events is 
of the order of 50 cm, which is the characteristic dimen- 
sion of the apparatus. That is, if these large-scale plumes 
originate from the boundary layer, it is as if the entire 
boundary layer on the bottom wall participates once in a 
while in this activity that we have called the large-scale 
plume. The maximum temperature in these large-scale 
plumes is a fraction of the excess temperature of the bot- 
tom plate (namely AT/2), so, presumably, the fluid that 
is participating in the formation of a typical large-scale 
plume comes from the top parts of the boundary layer; 
or, if a plume does indeed come from the very bottom 
parts of the boundary layer, it is already partially mixed 
by the time it reaches the probe midway between the 
top and bottom walls. They are certainly not small-scale 
events that scale on the thickness of the boundary layer. 
This description does not apply to small-scale plumes a-i, 
though, presumably, they too belong to the same family. 
In Fig. 1, we show the case of hot plumes, that is, the case 
when the wind at the measurement point arrives from the 
hotter bottom plate. One can imagine that the wind di- 
rection could be just the opposite, leading to the arrival 
at the probe of cold plumes coming from the colder top 
plate. We have analyzed such instances as well. Further, 
beyond a certain Rayleigh number, as described in Ref. 
U, the mean wind reverses itself randomly, so that a 
probe permanently held at one position sees hot plumes 
for some period of time and cold plumes for some other 
period of time. We have analyzed the two parts sepa- 
rately by stringing together only hot parts or the cold 
parts of a measured temperature trace. The data for the 
two cases have been examined separately. In each case, 
the telegraph approximation is related to specific proper- 
ties of the underlying physical processes associated with 
hot or cold plumes, as they encounter the probe during 
their motion. 
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FIG. 1: An example of the measured temperature time trace 
(upper part) and its random telegraph approximation (lower 
part). In addition to the small plumes marked a-i, there are 
several others such structures in the signals. They have not 
been marked merely because they occur below the zero line. 
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FIG. 2: Spectra of the telegraph approximation of two real- 
izations of the temperature trace compared with those of the 
temperature trace itself. The ordinate has been shifted so 
that the first points in both spectra coincide. 



II. RANDOM TELEGRAPH APPROXIMATION 

A more detailed discussion of the plumes will be pre- 
sented elsewhere but we limit ourselves here to a discus- 
sion of their tendency to cluster together occasionally. 
This is not obvious from the piece of the temperature 
trace shown in Fig. 1, and a longer trace crowds the 
plumes too much. It may be surmised that the clus- 
terization is indeed responsible for the mean wind in the 
apparatus. In the usual methods of analysis of turbulent 
signals ^4] , it is difficult to separate the clusterization ef- 
fect from the usual intermittency effects arising from am- 
plitude variability. To separate the two effects, we ignore 
the variation of the amplitude from one plume to another 
and replace the temperature trace of the type shown in 
the upper part of Fig. 1 by its random telegraph approx- 
imation, shown in the lower part. This approximation is 
generated from the measured temperature by setting the 
fluctuation magnitudes to 1 or depending on whether 
the actual magnitude exceeds the mean value (marked as 
zero and shown by the dashed line in the upper part of 
Fig. 1). Formally, for the temperature fluctuation 9(t) 
(with zero mean), the telegraph approximation T{t) is 
constructed as 




By definition, T can assume either 1 and 0. The tele- 
graph approximation can be generated by setting differ- 
ent "thresholds" than the mean. It turns out that most 



properties examined here are reasonably independent of 
the threshold; this comment will be made also at other 
specific places in the paper. It may be useful to know 
how the conventional statistics for the random telegraph 
approximation compare with those of the temperature 
signal. Figure 2 compares the spectral densities, E{f), 
of the telegraph approximations with those of the full 
signal. The comparisons are made separately for hot and 
cold cases. It is clear that the spectra of the telegraph 
approximation are close to those of the original signal in a 
significantly large interval of scales. The main difference 
is that the telegraph approximation is richer in spectral 
content above a certain frequency. This is not difhcult to 
understand from a visual inspection of Fig. 1. 

Of particular interest is the power-law behavior of the 
spectral densities of the telegraph approximation. For 
both hot and cold cases, they follow a power-law of the 
form 

E{f) - r". (2) 

Though this power-law behavior is clear from Fig. 2, we 
reproduce in Fig. 3 the spectra of the telegraph signal, 
computed from several records to attain better statistical 
convergence, in order to emphasize the power-law scaling. 
The exponent /? = 1.38 ± 0.02. A reasonable shift of the 
threshold does not change the spectral exponent (3. From 
the closeness of the spectra of the temperature trace with 
its telegraph approximation, it is inferred easily that the 
former has a spectral exponent of 1.38 as well, albeit over 
a smaller range of scales. Observation of the temperature 
trace spectra with such power law was first made in |^, 
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FIG. 3: Spectrum in the telegraph approximation computed 
using twenty realizations of the temperature trace, for the hot 
case (upper plot) and for the cold case (bottom part). The 
straight lines (the best fit) show that the power-law approxi- 
mation (2) holds well for both hot and cold cases. 
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FIG. 4: The probability density function of the duration r for 
the telegraph approximation of the temperature signal (for 
hot and cold cases). The straight lines (the best fits) are 
drawn to indicate the scaling law (3). 



and was explored theoretically in 6], and is now a well- 
known result. 

The probability density function of the duration be- 
tween events, r, for the telegraph approximation is shown 
in Fig. 4. Data for the hot and cold cases are given sep- 
arately. The log-log scale has been chosen to emphasize 
the power-law structure 



p{t) 



(3) 



For both hot and cold cases, we observe a = 1.37±0.03. 
A reasonable variation of the threshold in the vicinity of 
the average value does not change the exponent a. 

It is known that, for non- intermittent cases (Ref. 
and Sec. 4), the relation between the exponents a and /3 
is given by 



/? 3 - a. 



If we substitute in (4) the value of a ~ 1.37 (as observed 
in Fig. 4) we obtain /? ~ 1.63. This is considerably larger 
than that actual value measured in Fig. 3, namely 1.38. 
This discrepancy is the object of interest to us here; as al- 
ready remarked, since there are no amplitudes involved, 
it must be related to clusterization entirely (see also Q). 
It is of interest to note here that, for the telegraph ap- 
proximation of the temperature fluctuation in the tur- 
bulent atmospheric boundary layer, we have about the 
same values of a and /? as the present, while the spec- 
tral density of the temperature trace has a roll-off rate 
of about 1.66 (consistent with the Kolmogorov-Obukhov- 
Onsager-Corrsin theory ;4]). It is clear that the spectra 



of the temperature signal and its telegraph approxima- 
tion are closer in confined convection than in atmospheric 
turbulence. 



III. QUANTIFYING CLUSTERIZATION 

The difference between the observed telegraph spectral 
exponent {(3 ~ 1.38) and its value given by Eq. (4) (~ 
1.63) is a quantitative measure of clusterization (2j of 
plume-like objects observed in temperature traces. This 
is a part of intermittency. 

Intermittency of the so-called temperature dissipation 
rate 3,0 is characterized in turbulence by 



X 



1^1 



(5) 



(4) Following Obukhov 4], the local average 



Xt 



t + T 



X(t)dt 



can be used to describe the intermittency of x- The scal- 
ing of the moments. 



Ml 
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(6) 



assuming that scaling exists, is a common tool for the 
description of the intermittency Intermittent sig- 

nals possess a non-zero value of the exponent fiq. Of par- 
ticular interest is the exponent /i2 for the second-order 
moment. 
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This definition has a problem in the telegraph approx- 
imation because dT/dt is composed of delta functions. 
Fortunately, one is interested in the scaling of the dis- 
crete representation of the dissipation field, given by 

r 
n=l 

where n specifies an interval of space. This discrete defi- 
nition of Xt avoids the problem with delta functions. Ob- 
viously, for the telegraph signal, scaling exponents calcu- 
lated for the pulse-defmed dissipation, Eqs. (5)-(6), and 
for the discrete spike-like process are identical; this ob- 
servation is not true for the original temperature 
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FIG. 5: Normalized second moment of the local dissipation 
rate for the telegraph approximation plotted against r, for 
cold and hot cases. The straight lines (the best fits) are drawn 
to indicate the scaling law (6). 



The telegraph approximation is a composite of Heavi- 
side step functions, so the dissipation rate (5) is a com- 
posite of pulses (i.e. delta- functions) located at the edges 
of the boxes of the telegraph signal. For the uniform 
random distribution of the pulses along the time axis, 
fiq = 0. Non-zero values of /i^ mean that there is a clus- 
terization of pulses. Figure 5 shows dependence of the 
normalized dissipation rate (Xr) on r for the telegraph 
approximation of hot and cold signals. The straight lines 
(the best fits) are drawn to indicate the scaling law (6), 
for q — 2. It should be noted that scaling interval for the 
dissipation rate is the same as for the PDF and for the 
spectrum (cf. Figs. 3 and 4). Values of the intermittency 
exponent /12, calculated as slopes of the straight lines in 
Fig. 5 is /i2 ~ 0.47 ± 0.03. The relatively large value of 
the exponent ^2 suggests that the clusterization of the 
pulses is quite strong. The temperature dissipation can 
be also characterized by the "gradient" measure |jj,i8ji 



Xr 



where Vr is a subvolume with space-scale r (for detailed 
justification of this measure, see Ref. Q, p. 381 and 
later). The scaling law of the moments of this mea- 
sure are important characteristics of the dissipation field 
0. By Taylor's hypothesis we can replace dT/dx by 
dT/ (u)dt (where (u) is the mean wind and x is the co- 
ordinate along the direction of the wind), and can define 
the dissipation rate as 



Xr 



where r ~ r/{u). This, too, should follow the scaling 
relation (6). 



IV. CLUSTERIZATION AND THE SPECTRUM 

The spectrum of the telegraph approximation can be 
related to the probability distribution p(t) through 



Eif) = / W{t)Gi,ft)pit)dt, 



(7) 



where G{ft) is a transform function and the weight func- 
tion W{t) is supposed to have a scaling form 



Wit) - 



(8) 



Since the quantities G{ft) and p(t)dt are dimensionless, 
one can use dimensional considerations to find the expo- 
nent 6 (cf. Ref. II) through 



Wit) ^ (x) • t^ 



(9) 



and the use of (2), (3) and (7); this yields relation (4). 

To estimate the clusterization correction on the rela- 
tion between scaling spectrum and p(t), we should take 
into account two-point correlations in the telegraph sig- 
nal. This can be characterized by the two-point corre- 
lation function ^(i). In a situation where the two-point 
correlation function exhibits the scaling behavior 



(10) 



the correlation exponent 7 is the same as the exponent 
/i2. This is easily seen bv the well-known result that the 
correlation dimension Q D2 is related to 7 through 



^2 = 1-7, 



and to /i2 Ref. 2l "^^^ 



Do 



1 



(11) 



(12) 



thus yielding 7 = ^2- To estimate the weight function 
W{t) with the same dimensional considerations as above, 
and to take into account the two-point correlation (which 
characterizes clusterization), we replace (x) by 



t 



-M2/2 



(13) 



5 



Replacing (9) by 

mi)-(x?) (14) 

we have 

-t2-'^2/2. (15) 

The corresponding correction of Eq. (4) is 

/3 = (3 - /i2/2) - a. (16) 

Using the value of /i2 ~ 0.47 from Fig. 5 and the value 
a ~ 1.37 from Fig. 4, we obtain 

E{f) ^ (17) 

which compares well with the behavior found in Fig. 3 
(see also 0,0). 

V. SUMMARY AND DISCUSSION 

The results discussed so far are for a fixed Rayleigh 
number. The same situation seems to occur at other 
Rayleigh numbers. For those Rayleigh numbers where 
there is a reversal of the wind, the concatenated data 
corresponding to one given direction of the wind follow 
the same statistics as well. Thus, the characteristics 



discussed in the paper are generally valid for turbulent 
temperature fluctuations at all Rayleigh numbers cov- 
ered in the measurements. The main conclusion is that 
the telegraph approximation captures the main statis- 
tical features of the temperature time trace obtained 
in convection. This approximation, which gives a clear 
separation between clusterization and magnitude inter- 
mittency, has been useful in demonstrating that there is 
a significant tendency for the plumes to cluster together. 
The telegraph approximation turns out to be useful here 
because of the specific process of heat transport, which 
is determined in large measure by the random motion 
of temperature plumes. However, one can expect that 
this approximation (or its modifications) may be also 
useful in the description of other turbulent signals. The 
exponent /i2 for the telegraph approximation, completely 
determined by clusterization, is about 0.47. This should 
be compared with the intermittency exponent computed 
for the x~X correlation of the full temperature signal, 
which is about 0.36 — consistent with similar estimates 
available for passive scalars (see, e.g., Ref. 0). The 
clusterization exponent is thus larger than the classical 
intermittency exponent. From this, one can infer that 
the magnitude intermittency plays a smoothing role on 
the clusterization effects within the scaling interval. 

We thank L.J. Biven, J. Davoudi and V. Yakhot for 
useful discussions. 
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